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A fundamental prediction of our current hierarchical paradigm 
of galaxy formation is that massive galaxies form from an assem- 
bly of smaller units ^. The most massive objects form last, driven 
by the merging history of their dark matter haloes. The model 
succeeds in describing the clustering of galaxies^, but the evolu- 
tionary history of massive galaxies, as revealed by their visible 
stars and gas, presents problems. Near-infrared light allows us 
to measure the stellar masses of high-redshift galaxies^ and deep 
multi-colour images indicate that a large fraction of the stellar 
masses in massive galaxies form in the first 5 Gyr"*-^-^-^; but un- 
certainties remain due the lack of spectra to confirm the redshift 
and the role of obscuration. Here we report on the Gemini Deep 
Deep Survey, the deepest spectroscopic redshift survey ever un- 
dertaken which probes the most massive and quiescent galaxies 
back to an era only 3 Gyr after the Big Bang. We find that at 
least two thirds of massive galaxies have appeared since this era, 
but also that a significant fraction is already in place indicating a 
suprisingly slow rate of decline as we look back in time. 

Star- formation rates in the high redshift (z > 1) Universe have been probed 
by measuring rest-frame ultraviolet (UV) emission*^ of galaxies. The UV is dom- 
inated by newly formed massive stars but can be highly obscured^. Moreover 
star-formation in galaxies is a stochastic phenomenon and correlates poorly with 
mass: the most massive galaxies in the local universe are giant ellipticals and they 
have very weak UV emission 1° which makes them difficult to detect at z > 1. 
Direct determination of dynamical masses of galaxies requires spatially resolved 
velocity measurements of galaxies — this has been done out to 2 ~ ]^ii i2 ^^^^ 
such observations are very challenging at higher redshifts due to signal-to-noise 
limitations. 

An alternate approach is to infer the total mass in stars from light which traces 
accumulated stellar populations. Studies^ have shown that stellar mass correlates 
extremely well with dynamical mass out to z = 1. Stellar mass evolution can 
be predicted by using galaxy formation models based on cosmological numerical 
simulations augmented with analytical star-formation recipes,-- which adopt sim- 
plified prescriptions for various heating and cooling processes in the interstellar 
mediumi^. 



2 



Deep near-infrared data provides a window on stellar masses of galaxies at 
high redshift. This light is dominated by the old, evolved stellar populations in 
galaxies and is little affected by transient star-formation^^ or dust obscuration. 
To a good approximation the total i^'-band light traces the accumulation of stel- 
lar mass; equivalently one can say that the stellar mass-to-light ratio (M/Lk) 
is nearly constant with little dependence on the previous star-formation history 
(SFH). In fact rest- frame M/ Lx varies only by a factor of two between extremely 
young and extremely old galaxy stellar populations, in contrast M/ Lb can vary 
by more than a factor of ten and M/ Luv varies tremendously (reaching as low 
as 1% of the solar value) because the UV light is dominated by the instantaneous 
star-formation rate (SFR) rather than the stellar mass. Even at redshifts up to 
z = 2 the observed-ixdimQ K-band probes the rest-frame i?-band which is still 
dominated by old stellar populations. The M/Lr variation is still less than a 
factor of four, partly due to the youth of the Universe at that redshift. 

We use the Gemini Deep Deep Survey (GDDS)^^ to define a new sample of 
150 galaxies with K < 20.6 and 0.8 < z < 2 located in four independent 30 
arcmin^ fields. The spectroscopic identification completeness is 89%. The de- 
termination of stellar masses for each galaxy follows standard multi-color stellar 
population fitting techniques; this approach is fairly general and robust and is de- 
tailed in the Methods appendix. We find that the i^-band light traces the stellar 
mass quite well; evolutionary changes in galaxy numbers with redshift are much 
more important than changes in color. The findings presented below are simply 
driven by the presence of numerous ~ 20 galaxies at 2; > 1.5 which must be 
massive objects. 

Figures 1 & 2 illustrate the nature of the most massive galaxies. They are 
found by the ODDS to 2; = 2 at ~ 20, even when they are very red (I — K > 
4). The red galaxies are predominantly red due to old stellar populations (these 
are GDDS spectral classes^'' '001' and 'Oil' showing photospheric features from 
evolved stars) and not dust reddening of young lower mass galaxies. The red, old 
galaxies make a large contribution (30%) to the stellar mass in the Universe in the 
redshift range 1.2 < z < 1.8. Also at 2; > 1 there are a number of blue galaxies 
despite this being a i^-selected sample. These correspond to massive star-forming 
galaxies which have high metal abundances 

The cumulative stellar mass density per unit volume in each redshift bin down 
to various mass thresholds (Figure 3; Table 1) is computed following the stan- 
dard V/Vmax formalism (this corrects for the smaller redshift ranges covered by 
fainter objects ^^) for a K < 20.6 limit and weighting by the sampling. We use 
K-corrections from our individual spectral energy distribution (SED) fits, but the 
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results do not depend strongly on the details of the K-correction. Although a K- 
selected sample is a good proxy for a mass- selected sample, it is still necessary to 
consider incompleteness as a function of M/L. To do this we compute the max- 
imum possible M/L at each redshift for a model galaxy as old as the Universe 
which formed all its stars at once (a 'Simple Stellar Population' or SSP). This 
M/ L is converted to a mass limit (via the i^'-band flux limit and the K-correction) 
and is shown in Figure 1. We are complete above this mass limit; bluer objects 
can be seen below this limit. Bins which might miss high M/L objects are plot- 
ted as lower limits in Figure 3. The error bars are calculated from shot noise 
on the number of galaxies in each bin. Of course these do not include the ef- 
fects of large-scale structure, but we believe these are not significant because our 
GDDS fields are large, were selected from even larger area images in regions near 
average density'^ (i.e. neither highly over-dense nor under-dense) and because 
our colour-dependent weights normalise to the full imaging area (554.7 arcmin^) 
which would counteract the additional clustering of red objects. Finally we have 
performed the check of splitting the sample by different fields (different indepen- 
dent sight-lines). The same general results are found for these albeit with larger 
errors. We have also assessed the effect of the mass fitting errors on the mass den- 
sities with our Monte-Carlo methods; this is not a significant source of error. In 
every bin galaxies with spectroscopic redshifts dominate the mass budget except 
for the 1.6 < 2 < 2 and M > 10^^-^Mq bin, where the spectroscopic complete- 
ness is only 50% and we augment our spectroscopy with photometric redshifts. 
Analysing the sub-sample with only spectroscopic redshifts results in no signif- 
icant changes to any bin except for this one (which is thus 0.3 dex lower); the 
essential scientific result is unchanged. 

There is a clear decline in stellar mass locked up in the most massive galax- 
ies. This trend is in accord with ideas of gradual assembly. However it de- 
clines only slowly towards high redshift and, surprisingly, the massive galaxies 
do not decline more rapidly than the whole population. We note that galaxies with 
M > 10^°-^ Mq, are brighter than our flux limits (for any possible M/L value) 
throughout our redshift range; a regression on these gives an acceptable fit for 
p oc (1 + z)-^-'^'^'^-^. At ;z = 1 the mass densities for the M > IO^^-^Mq sample 
are 38 (±18)% of their local value^"; at z = 1.8 this becomes 16 (±6)%. These 
results for the most massive galaxies are consistent with previous Hubble Deep 
Field South photometric redshift determinations"^-^, but inconsistent with the more 
rapid decline (factor of ~ 6 over < z < 1) found with photometric redshifts 
by the Munich Near- Infrared Cluster Survey^ (however we note that their SED 
fitting forces maximally old galaxies at all redshifts, which could be problematic). 
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Overlaid on Figure 3 (shaded region) is the range of estimates for the growth of 
total stellar mass based on the integral of the observed rest-UV derived SFR-^ 
relationship. (These are based on the an analytic fit^° of points from Figure 9 of 
Steidel et al.^, both with and without extinction correction, integrated using PE- 
GASE.2^^) Our uppermost points represent lower limits on the total stellar mass 
density and we confirm previous findings that an extinction correction is essential 
for UV SFR estimates to be consistent with stellar mass measurements^. 

Theoretical models of galaxy formation convert gas into stars in dark matter 
halos using semi-empirical recipes and must satisfy three key observables: the 
first is the distribution of local galaxy luminosities, the second is that of galaxy 
colors and the third is the abundance of massive galaxies at high redshift. The last 
of these has until now been the most difficult to measure from observations. In 
Figure 3 we plot the abundance predictions of the 'GALFORM' models '"^-^^ which 
satisfy the first two constraints; it is evident in this model that massive galaxies 
disappear much more rapidly than we see in our data; this is because the model 
stellar mass build-up traces the merging of cold dark matter haloes. In particular 
GALFORM has a strong dependence of evolutionary rate on mass which is not 
seen in our data. A similar theoretical comparison was made on the Hubble Deep 
South photometric data with similar conclusions^. 

We note that our measured abundance of massive galaxies at 2; ~ 2 does not 
violate the boundary constraints of cold dark matter models, by which we mean 
taking the predicted abundance of 2; = 2 massive dark matter halos and scaling 
by the cosmological baryon / dark matter density^ (i^b/fi^ = 0.17). This gives 
a predicted baryonic mass density of 10^ Mq Mpc"'^ in haloes with baryonic 
masses > lO^^M©, a factor of ten above our measurements. Thus the models can 
match our densities if only ~10% of the baryons in massive haloes are converted 
in to stars hy z = 2. We note though that in the Universe today, the stellar den- 
sity ^° = 0.004 gives r2*/i7fe = 0.1; this would imply massive haloes at z = 2 
have managed to convert baryons into stars with the same overall efficiency as the 
average Universe had achieved by 2 = 0. Massive haloes must have much greater 
star- formation efficiencies at earlier times. 

Models have been proposed which adjust the star-formation histories in 
this way and which would be more in accord with our findings of a high abun- 
dance of massive galaxies at early times. However they fail to match existing data 
on galaxy colours ^^-^^ or the galaxy luminosity function— so we defer detailed 
comparison to future work. 

What is observationally clear is that we have measured directly the abundance 
of massive galaxies out to z = 2 for the first time from a deep and relatively wide- 
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area survey with secure redshifts and spectroscopic classifications for galaxies. 
We find (with 99% confidence) that at least two-thirds of the mass in these objects 
has formed since z = 1.8. However we do not find a more rapid evolution of 
the giant population compared to the Universal average as would be expected if 
stellar mass grew simply proportional to dark matter assembly. At 2; < 2 we 
find much of the mass is in galaxies with old stellar populations; these objects 
are plausible precursors of modern massive elliptical galaxies. This new data is 
broadly consistent with developing ideas that models of galaxy formation have 
to be tuned towards earlier star-formation in high-mass haloes. The abundance 
of massive galaxies at z < 2 is now firmly established and will provide the vital 
missing leg in the tripod of key observations required to understand how galaxies 
form. 

Methods 

Stellar Mass Fitting 

We adopt a cosmology^ of fi^ = 0.3, f^A = 0.7, Hq = 70 kms~^Mpc^^ and 
use Vega magnitudes. For the mass-function analysis we use GDDS galaxies with 
K < 20.6 and 0.8 < z < 2 (150 galaxies; 89% spectroscopic completeness based 
on Abraham et al. — identification confidence classes >= 2; galaxies without 
spectroscopic redshifts are assigned to their photometric redshifts). To account 
for the higher priority given to red galaxies in the slit-mask design we use the 
'sampling weights' from Abraham et al. which are the selection probability for 
spectroscopy as a function of I — K and K. 

For each galaxy we derive the most likely stellar mass and a range of uncer- 
tainty. These are derived by evaluating against model optical-IR spectral-energy 
distributions (SED) to determine the M/L^ ratio, and hence the mass. We note 
that a variety of approaches to accomplish this have been described in the literature 
which vary in the level of detail in which they treat star formation. For example. 
Cole et al. ^° used a simple set of monotonic star- formation histories (SFH) with 
a varying e-folding time-scale and a fixed dust law. A potential problem is posed 
by the fact that real galaxies have more complex SFHs, for example a recent star- 
burst can make an old galaxy temporarily bluer and lead to an underestimate of 
M/L using this method. One approach to account for this effect is to introduce a 
second young SED component^ superimposed on the old population; this can be 
computationally expensive depending on the amount of freedom allowed for the 
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second component. 

In a spectroscopic sample there is no degeneracy between SED model fitting 
and photometric redshift (which is also based on SED fitting). We adopt two- 
component modelling (using PEGASE.2^^ to calculate spectra) in order to be able 
to assess the possible biases due to starbursts on the calculated masses and al- 
low a range of dust extinction (0 < Ay < 2 mag) and metallicity (0.0004 < 
Z < 0.02). The primary component is modelled using a star-formation rate 
SFR oc exp(-t/r) with r =0.1, 0.2, 0.5, 1, 2, 4, 8 & 500 Gyr (the first ap- 
proximates an instantaneous starburst and the last a constant SFR). The secondary 
component is a starburst, modelled with a r = 0.1 Gyr exponential, which can 
occur at any time and have a mass between 10 ~^ and twice that of the primary 
component. The model age is constrained to be less than that of the Universe at 
the appropriate redshift, but any formation epoch is allowed. Like all studies of 
high redshift star-formation, we must assume a Universal Initial Mass Function 
(IMF). We primarily use the BG03 IMF^' which fits local cosmic luminosity den- 
sities well; it has a similar high-mass slope as the classical Salpeter IMF^^ but 
with a more realistic break at O.5M0. Stars with masses below the break never 
contribute significantly to optical/IR light so the overall effect is to re-scale the 
total stellar masses to more reasonable values. Similarly we re- scale literature 
Salpeter based numbers '^■^■2" to the BG03 IMF using Msp = 1.82Mbg- Our re- 
sults are robust for any reasonable choice of IMF with a similar slope; to illustrate 
this we have also calculated masses using the popular Kennicutt IMF—. 

Our approach is to find via exhaustive grid search all models consistent with 
the Viz' K photometry in the observed frame of each galaxy. This colour set is 
available for all galaxies and covers rest-frame ultraviolet through near-infrared. 
We do not include the actual spectra (apart from the redshift information) in the 
fits because of variable quality and signal-to-noise, however we find the spectral 
classes are broadly consistent with the best photometric SED fits. The full distri- 
bution function of allowed masses were calculated by Monte-Carlo re-sampling 
the photometric errors. The final masses and error bars represent the mean and 
standard deviation of this full distribution function. Typically, we find the masses 
are fitted to ±0.17 dex in the K < 20.6 sample. The stellar masses are very robust 
against the details of the fitting. Using the Monte-Carlo machinery to investigate 
the effect of different assumptions about metallicity, dust and bursts we find that 
the largest effect is due to bursts. If we disallowed bursts then the masses typi- 
cally decrease by only 0.2 dex. Finally we note the variation in M/ Lk over the 
range 1 < 2; < 2 is constrained by the age of the Universe (6-3 Gyr), typically the 
maximum range in our sample is a factor of three. 
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Table 1: Stellar Mass Density measurements in the universe as a function of 
galaxy mass threshold and redshift. p is the stellar mass density in galaxies with 
stellar masses > Mum in Mq Mpc"^. The values of phi and pio are the la upper 
and lower limits from counting statistics. BG03 IMF mass thresholds and mass 
densities can be multiplied by 1.82 to convert to the Salpeter IMF. 
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Figure 1: Stellar mass-redshift distribution for our galaxies. The 1cr error 
bars come from the IVIonte-Carlo mass fitting. Symbol colours code the ob- 
served K-magnitude (see panel). Solid/open symbol shapes denote spec- 
troscopic/photometric redshifts respectively. Circle/star symbols denote 
objects redder/bluer \nV - I than a model Sbc galaxy template ^° respec- 
tively. Symbol size is l^eyed to the /-band magnitude. The horizontal line 
denotes the characteristic Schechter mass scale in the local Universe^°.. 
The solid curve shows how the K-i\ux limit translates into a mass com- 
pleteness limit for a maximally old Simple Stellar Population. The dashed 
curve shows an example mass limit for bluer objects (SFR=const. model). 
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Figure 2: Colour-redshift distribution for our galaxies. Observed frame 
I - K colour is plotted with symbol size keyed to mass, shape to op- 
tical V - I colour (as in Figure 1) and symbol colour l<eyed to the 
spectral classification (see panel). Solid/open symbols denote spectro- 
scopic/photometric redshifts. IVIodel tracl<s (solid lines) are shown rang- 
ing from a a maximally old Simple Stellar Population (SSP) through to 
SFR=const., for synthetic solar metallictiy galaxies which form at 2; = 10. 
(Note the points red-ward of the old tracl< at 2; = 0.6 may be artificially red- 
dened by observational aperture/metallicity effects, this is insignificant for 
z > 0.8.) 



05 - 



n 

a 

o 

1^ 



A 



00 

o 



CD - 



/SFR (raw & corrected) 




• GDDS 


logio(M) > 10.2 


oCole et al. 2001 


logio(M) > 10.5 


GALFORM 


logio(M) > 10.8 


□ Dickinson et al. (2003) 


TOTAL 


-Fontanaet al. (2003) 


9.7<logi(,(M)<10.7 


1 . , , . 1 , 


0.5 


1 




Redshift 







1.5 



Figure 3: Mass density in stars vs redshift. Values from our spectroscopic 
sample are compared with previous estimates which are all based on pho- 
tometric redshifts"^-^ (except for the local z = 0.1 point^°). We plot the 
cumulative mass density of galaxies (converted to our Initial Mass Func- 
tion choice; error bars are la) more massive than a given mass threshold 
(see legend for threshold colour coding). Theoretical 'GALFORM' models 
are also plotted (see text). The shaded region shows the result of integrat- 
ing the Universal ultraviolet-derived star-formation history with and without 
a dust extinction correction. 



